Origin of Spin Ice Behavior in Ising Pyrochlore Magnets with Long 
Range Dipole Interactions: an Insight from Mean-Field Theory 
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Recent experiments suggest that the Ising pyrochlore magnets Ho2Ti2C>7 and Dy2Ti2C>7 display 
qualitative properties of the ferromagnetic nearest neighbor spin ice model proposed by Harris et 
al, Phys. Rev. Lett. 79, 2554 (1997). The manifestation of spin ice behavior in these systems de- 
spite the energetic constraints introduced by the strength and the long range nature of dipole-dipole 
interactions, remains difficult to understand. We report here results from a mean field analysis that 
shed some light on the origin of spin ice behavior in (111) Ising pyrochlores. Specifically, we find 
that there exist a large frustrating effect of the dipolar interactions beyond the nearest neighbor, and 
that the degeneracy established by effective ferromagnetic nearest neighbor interactions is only very 
weakly lifted by the long range interactions. Such behavior only appears beyond a cut-off distance 
corresponding to O(10 2 ) nearest neighbor. Our mean field analysis shows that truncation of dipolar 
interactions leads to spurious ordering phenomena that change with the truncation cut-off distance. 
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I. INTRODUCTION 

The past five years have seen a resurgence of signif- 
icant interest devoted to the systematic study of geo- 
metrically frustrated magnetic systems . Frustration 
arises when a magnetic system cannot minimize its total 
classical ground-state energy by minimizing the energy of 
each spin-spin interaction individually 0. This most of- 
ten occurs in materials containing antifcrromagnetically 
coupled magnetic moments that reside on geometrical 
units, such as triangles and tetrahedra, which inhibit the 
formation of a collinear magnetically-ordered state. 

It is very common for models of highly frustrated mag- 
netic systems to display ground state degeneracies where 
the system is "underconstrained" , giving many distinct 
spin configurations give the same ground state energy. 
The best known example is probably that of the two- 
dimensional Ising antiferromagnet on a simple triangu- 
lar lattice, for which it was shown by Wannier that the 
ground state is macroscopically degenerate with an ex- 
tensive ground state entropy . Another example is the 
nearest-neighbor antiferromagnet face-centered cubic lat- 
tice with Heisenberg [p|-|l2|], and Ising spins (jl3|. In the 
past ten years, much attention has been devoted to the 
kagome lattice of corner-sharing triangles, and the py- 
rochlore lattice of corner-sharing tetrahedra with Heisen- 
berg spins interacting via nearest neigbor antiferromag- 
netic interactions (j]-^]. 

Typically, there are two generic class of mechanisms 
that can lift the degeneracies in frustrated magnetic sys- 
tems. The most common and simplest mechanism pro- 
ceeds via energetic perturbations beyond the term(s) in 
the spin Hamiltonian that cause the degeneracies. De- 
generacies can also be lifted by either thermal fluctu- 
ations or quantum zero-point fluctuations |j|||,[l4]dl9]] . 



There are a number of well known situations where spin 
interactions beyond nearest neighbor cause long range or- 
der that would otherwise be absent in a highly frustrated 
lattice with only nearest neighbor antifcrromagnetic ex- 
change present. For example, the classical kagome p(J 
and pyrochlore Jl6| , p~7 21 nearest neighbor Heisenberg 
antiferromagnet lattices display either Neel order or a 
dramatic reduction of ground state degeneracy when ex- 
change interactions beyond nearest neighbor are con- 
sidered. Long range dipole-dipole interactions in the 
FCC Heisenberg antiferromagnet Jll],|l2| and in the py- 
rochlore lattice with Heisenberg spins (as pertains to the 
Gd2Ti207 pyrochlore material |2g-gJ|), have also been 
shown to reduce significantly the degeneracy otherwise 
present in the nearest neighbor version of these systems. 

While high geometric frustration usually arises in anti- 
ferromagnetic systems, Harris and collaborators recently 
showed that the pyrochlore lattice of corner sharing tetra- 
hedra with Ising spins pointing along a local cubic (111) 
axis (See Fig. 1) constitutes an interesting and unusual 
example of high geometric frustration when the near- 
est neighbor interaction is actually ferromagnetic |^5|,^6| . 
Harris and collaborators introduced the concept of spin 
ice to emphasize the analogy between ferromagnetically 
coupled (111) Ising moments on the pyrochlore lattice 
and the problem of proton ordering in common hexagonal 
ice, Ih ]25|-|2^|. In the simple model of nearest neighbor 
ferromagnetic interactions, the system has the same 'ice 
rules' for the construction of its gro und state as those 
for the ground state of real ice p9), hence the name 
spin ice. The resulting extensive ground state entropy of 
hexagonal ice was first estimated by Pauling in 1935 ]3C| ] 
which, remarkably, provides a fairly accurate estimate of 
the residual low temperature entropy of real pyrochlore 
spin ice materials |}1 



1 




FIG. 1. The lower left 'downward' tetrahedron of the py- 
rochlore lattice shows Ising spins (arrows). Each spin axis 
is along the local (111) quantization axis, which goes from 
one site to the middle of the opposing triangular face (as 
shown by the disks on the triangular faces) and meets with 
the three other (111) axes in the middle of the tetrahedron. 
For clarity, black and white circles on the lattice points denote 
other spins. White represents a spin pointing into a downward 
tetrahedron while black is the opposite. The entire lattice is 
shown in an ice-rules state (two black and two white sites 
for every tetrahedron). The hexagon (thick gray line) shows 
a low energy "loop" excitation (see text) which corresponds 
to reversing all colors (spins) on the loop to produce a new 
ice-rules state. 

Interestingly, materials with local (111) Ising spins 
are realized in Ho 2 Ti 2 7 [|§(3|^|, Dy 2 Ti 2 7 |l| and 
Tb 2 Ti 2 07 |35|-|3^|, where the rare earth magnetic mo- 
ments (Ho^7 Dy 3+ and Tb 3+ ) reside on the sites of 
the pyrochlore lattice. For (111) Ising pyrochlore sys- 
tems, each moment points along the axis joining the 
centers of the two tetrahedra that it belongs to (see 
Fig. 1). For ferromagnetic nearest neighbor interaction, 
the ground state is macroscopically degenerate, but with 
the constraint that two moments must point in and two 
must point out of every tetrahedron which, as mentioned 
above, is a constraint that maps precisely onto the ice 
rules j25]-§0). 

The nearest neighbor model of Harris and Bramwell 
p6[ shows no evidence of a magnetic transition to long 
range magnetic order |p7| , p8| , ^9| , but displays a broad 
peak in the specific heat at a temperature of the or- 
der of the nearest neighbor ferromagnetic coupling. Both 
Ho 2 Ti 2 7 H,||-|il] and D y2 Ti 20 7 |l| show qualitative 
behavior consistent with this spin ice picture. 

As explained above, the occurrence of a macroscopic 
degeneracy in frustrated magnetic is due to "under- 
constraints" ; many different spin configurations (here 
the ice-rules "two-in— two-out" per tetrahedron) give the 
same ground state energy We therefore expect that 
the microscopic origin of the spin ice phenomenon in 
Ho 2 Ti 2 7 and Dy 2 Ti 2 7 has to directly involve such a 
manifestation of underconstraints in the energetics in- 
volved in the formation of the spin ice ground state. 



Our goal is therefore to understand how the competing 
energy scales in the spin ice materials are able to produce 
underconstraints. In order to identify such undercon- 
straints, we must first consider the various energy scales 
at play in the above materials. Firstly, the single-ion 
ground state for Ho 3+ and Dy 3+ in the pyrochlore struc- 
ture is well described by an effective classical Ising dou- 
blet with a ground state magnetic moment of /i « 10/iB 
and a nearest neighbor distance, r nn , of approximately 
3.54A for both materials [[33). For two nearest neighbor 
moments pointing along their local (111) direction, the 
nearest neighbor dipole-dipole energy scale, D nn , is 

D m =\{^)4- «+2-35K , (1.1) 

where, as discussed in the next section, the 5/3 fac- 
tor comes from the orientation of the Ising quantization 
axes relative to the vector direction connecting interact- 
ing nearest neighbor magnetic moments. Interestingly, 
D nn is positive (i.e. ferromagnetic) and therefore, as in 
the Harris and Bramwell nearest neighbor model |2q] , is 
frustrating for Ising spins on a tetrahedron. 

The experimentally determined Curie- Weis tempera- 
ture, $cWj extrapolated from temperatures below T ~ 
100K is +1.9 K for Ho 2 Ti 2 7 || and +0.5 K for 
Dy2Ti 2 7 fnl , respectively. These two values show that 
$cw is of the same order of magnitude as the near- 
est neighbor dipolar energy scale. Furthermore, it is 
well known that rare-earth ions possess very small ex- 
change energies. Consequently, dipole-dipole interactions 
in Ho 2 Ti 2 7 and Dy 2 Ti 2 7 constitute a preponderant 
interaction, as opposed to a weak perturbation as in mag- 
netic transition elements where the exchange interaction 
predominates over dipolar interactions. As we discuss 
below, fits to experimental data show that the nearest 
neighbor exchange interaction is antifcrromagnetic for 
both Ho 2 Ti 2 7 and for Dy 2 Ti 2 7 and, when considered 
alone, is therefore not frustrating for (111) Ising moments 
on the pyrochlore lattice fl^J^Jlof . In order to consider 
the combined role of exchange and dipole-dipole inter- 
actions, we find it useful to define an effective nearest 
neighbor energy scale, J e s, for (111) Ising spins: 

Jeff = Jnn + Dnn , (1-2) 

where J nn is the nearest neighbor exchange energy be- 
tween (111) Ising moments. This simple description pre- 
dicts that a (111) Ising system would display spin ice 
properties even for antifcrromagnetic nearest neighbor 
exchange, J nn < 0, but has long as J c ff = J nn + -Dnn > 
0. Fits to experimental data give J nn ~ —0.52 K 

for Ho 2 Ti 2 7 @ and J nn 1.24 K for Dy 2 Ti 2 7 

|fll| . Thus, J c ff is positive (using D nn — 2.35K) hence 
ferromagnetic and frustrated for both Ho 2 Ti 2 7 and 
Dy 2 Ti 2 7 , but not for Tb 2 Ti 2 7 for which J nn « -0.9 
K and D nn = +0.8 K |35| [38|. It would therefore appear 
natural to ascribe the spin ice behavior in both Ho 2 Ti 2 7 
and Dy 2 Ti 2 7 to the positive J c g value as in the simple 
model of Harris and Bramwell |^,^6|. However, the sit- 
uation is much more complex than it naively appears, 
especially when considering the spin ice phenomenon in 
the context of an underconstrained problem. 
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The main conceptual difficulty in intuitively under- 
standing the physical origin of spin ice behavior in 
rare-earth titanates stems from the very nature of the 
large dipole-dipole interactions in the Dy2Ti207 and 
Ho2Ti207 materials. Recall the general discussion above 
on the role of perturbations in frustrated magnetic sys- 
tems with degenerate ground states. Dipole-dipole inter- 
actions are "complicated" in that (i) they are strongly 
anisotropic since they couple the spin, S*% and space, 
Yij, directions, and (ii) they are also very long range 
(oc 1 /rfj ) . One would naively expect that both the spin- 
space coupling, (S^ 1 • Vij), and the long range nature 
of the interaction (beyond nearest neighbor) introduces 
"so many" energetic constraints on the ground state spin 
correlations that consideration of the dipolar interaction 
beyond nearest neighbor would lift most, if not all lo- 
cal degeneracies. Such an effect would give rise to a 
long range ordered state at a well defined (possibly in- 
commensurate) wave vector at some critical temperature 
T C (J,D) ~ 0(D nn /kB) ~ 2 K, as opposed to spin ice 
behavior. This simple argument would appear to be fur- 
ther supported by noting that recent calculations show 
that vanishingly small, but nonzero long range dipole- 
dipole interactions select a unique (non-degenerate) long 
range Neel ordered state for Heisenberg spins in an other- 
wise classical nearest neighbor pyrochlore antiferromag- 
net p2| , p3[ . This clearly does not happen in real spin ice 
materials such as Ho 2 Ti 2 7 HQ and Dy 2 Ti 2 7 fl), 
at least down to a temperature T ~ 200 mK. Hence the 
question that we address here is: 

When there is an effective ferromagnetic nearest neighbor 
interaction, J Q g, in (111) Ising pyrochlores, why do long 
range dipolar interactions fail to destroy spin ice behav- 
ior, and not, instead, give rise to long range Neel order 
with a critical temperature T c ~ 0(D nn )? 

A simple explanation to this question might be that 
the inability of long range dipolar interactions to lift 
the ground state degeneracy established at the nearest 
neighbor level arises from a "mutual frustration" of the 
degeneracy-lifting energetics, or mean-field, coming from 
the spins beyond the nearest neighbor distance. This ar- 
gumentation, if correct, suggests that a mean- field theory 
which contains dipolar interactions up to an "appropriate 
cut off" distance, and aimed at determining the ordering 
wave vector as a function of the cut-off distance of the 
dipolar interactions may be used as a first step to investi- 
gate the microscopic origin of spin ice behavior, and the 
failure of a dipolar-driven degeneracy lifting process at a 
temperature T ~ 0(D nn ). As we show below, we find 
that the truncation of dipolar interactions lead to spu- 
rious results when the truncation distance is less than 
< 100 nearest neighbors, and that spin ice behavior is 
restored only when considering the dipolar interactions 
to very large cut-off distances. 

The rest of the paper is organized as follows. In the 
next section we present a description of the mean-field 
theory we use to determine the ordering wavevector in 
(111) Ising pyrochlores as a function of cut-off distance. 
Our results are presented in Section III, followed by a 



brief discussion in Section IV. 



II. MEAN-FIELD THEORY 

Our aim in this section is to use mean-field theory 
to determine the the critical ("soft") modes and, con- 
sequently, the nature of the magnetic phase(s) exhibited 
by a classical model of Ising spins on a pyrochlore lattice 
with local axes (the (111) directions of the cubic unit 
cell). This system is described by a Hamiltonian with 
nearest neighbor exchange and long range dipolar inter- 
actions f[4ll|42[: 



(ij) 



Dr A 



E 



Sf 4 • s* 



3(81 



(2.1) 



The first term is the near neighbor exchange interaction, 
and the second term is the dipolar coupling between the 
(111) Ising magnetic moments. For the open pyrochlore 
lattice structure, we expect very small second and fur- 
ther nearest neighbor exchange coupling Ej. We there- 
fore only consider exchange interactions between near- 
est neighbor spins. Here the spin vector S*' labels the 
Ising moment of magnitude \S*'\ — 1 at lattice site i 
and oriented along the local Ising (111) axis The 
distance |ry| is measured in units of the nearest neigh- 
bor distance, r nn . Here J represents the exchange en- 
ergy and D = (/i /47r)/Lt 2 /rjJ n . Because of the local 
Ising axes, the effective nearest neighbor energy scales 
are J nn = J/3 and, as mentioned above, D nn = 5-D/3, 
since Zi ■ Zj = —1/3 and (Bi ■ ry)(ry • Zj) — —2/3. 

We now proceed along the lines of Reimers, Berlin- 
skly and Shi in their mean-field study of Heisenberg py- 
rochlore antiferromagnets plj ], and which has recently 
been extended to investigate systems with long range 
dipole-dipole interactions [^2|,^3) . We consider the mean- 
field order parameters, B(ri) at site r^. The pyrochlore 
lattice is a non-Bravais lattice, and we use a rhombohc- 
dral basis where there are four atoms per unit cell located 
at (0, 0, 0), (1/4, 1/4, 0), (1/4, 0, 1/4), and (0, 1/4, 1/4) in 
units of the conventional cubic unit cell of size a = r nn <\/8. 
Each of these four points define an FCC sublattice of 
cubic unit cell of size a. We relabel the spins, S(rj), 
in terms of unit cell coordinates, and a sublattice index 
within the unit cell, and take advantage of the trans- 
lational symmetry of the lattice by expanding the order 
parameters B(ri) in terms of Fourier components. In this 
case B a (r;) = B(ri)zf on the a'th sublattice site of the 
unit cell located at r, can be written as 



S°(r i )5? = y;B B (q)ifexp(iq.r i ) 



q 



(2.2) 



where zf is a unit vector along the local (111) Ising axis 
on the i'th site of sublattice a. The spin-spin interaction 
matrix, l 7 ab (|r i j |), including both exchange and dipolar 
interactions, reads: 
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J a \\v ij \) = J{zt-z))8 rij 



(2.3) 



D, 



dd ■ 



2a j.ab ~b j.ab | 

■■>> ' " ^ , (2.4) 



where 8 a a is the Kronecker delta, rff denotes the inter- 



spin vector Yij that connects spin S° on the a sublattice 
ltfj 

of its Fourier components as 



to spin S b on the b sublattice. We write J ab (\Yij | in terms 



J 



ab ( 



(2.5) 



where M ce u is the number of unit cells with 4 spins per 
unit cell. The quadratic part of the mean-field free- 
energy, F^ 2 \ then becomes 21 1: 



F^(T)/M cell = \ J> a (q) {T6 ab - J ab (q)} B b (-q) 



q,(a&) 



(2.6) 



where T is the temperature in units of 1/fce- Diagonal- 
izing F^ (T) requires transforming to normal modes of 
the system 



B a (q) 



(2.7) 



where {$ Q (q)} are the eigenmodes, and t/(q) is the uni- 
tary matrix that diagonalizes J ah ((\) in the sublattice 
space, with eigenvalues A Q (q) 



J2j ab (q)U ba (q) = \ a (c l )U aa (q) 



(2.8) 



Henceforth we use the convention that indices (ab) label 
sublattices, and index a labels the normal modes. We 
express F^ 2 \T) in terms of normal modes as 

F^/M cell = 1^5>"(q)<I> Q (-q){fcBT- X a (q)} . 



(2.9) 

An ordered state first occurs at the temperature 

T c = max q {A mQa; (q)} , (2.10) 

where X max (q) is the largest of the four eigenvalues 
(a = 1,2,3,4) at wavevector q, and where max q indi- 
cates a global maximum of the spectrum of X max (q) for 
all q. The value of q for which A" (q) is maximum is the 
ordering wavevector q or d- 

Let us briefly explain how we proceed using the above 
set of equations to determine the critical ( "soft" ) mode(s) 
of the system at T c . The Fourier transform of J ab (|rjj|) 
is calculated using the inverse transform of Eq. (2.5). 
The pyrochlore lattice has a symmetry of inversion with 
respect to a lattice point and this implies that J ab (q) 
is real and symmetric. The eigenvalues and eigenvectors 
are found using a standard numerical packages for eigen 
problems of real symmetric matrices. 



III. RESULTS 

For each q there are 4 eig envalues of J ab {q) (X a (q) , 
a = 1, 2, 3, 4). To determine q or d we need to find the 
q value for which A" is maximum. We therefore calcu- 
late X max (q), the largest of the 4 eigenvalues of J ah (q) 
for each q. Figures 2—6 show X max (q) vs q in the (hhl) 
plane for various cut-off distance, N c , of the dipolar inter- 
actions used in t he c alculation of J ab (q) via the inverse 
transform of Eq. 2.5. Figures 2—6 correspond to cut-off 
distances for the 1st, 2nd, 5th, 100th and 1000th near- 
est neighbor (N c = 1, 2, 5, 100 and 1000), and which 
correspond to physical distances R C (N C = l)/r nn = 1, 
R C (N C = 5)/r nn = y/7, R C (N C = 10)/r nn - 
R C (N C = 100)/r nn = VI36 and R C {N C = 1000)/r nn = 37, 
expressed in units of the nearest neighbor distance, r nn . 
In the calculations here, we have used J = and measure 
X max (q) in units of the dipolar strength, D. Our conclu- 
sions below are independent of the choice made for J as 
long as we are in the spin- ice regime (J/3 + 5D/3) > 0. 
For J sufficiently negative (antiferromagnetic) the or- 
dering wavevector is at q or d = 0, corresponding to a 
two-fold Ising state all-in— all-out as discussed in Refs. 
|l8| , p6 40, 4l}|. Several comments are in order. 

Cutting off the dipolar interaction at N c — 1 gives the 
effective nearest neighbor Hamiltonian: 



H eS (N c - 1) 



5_D 

IT 



(3.1) 



where a. t = ±1. For J cS = (J + 5L>)/3 = J nn + D nn > 0, 
and we recover the nearest neighbor Ising spin ice model 
of Harris and Bramwell p5| , p6[ . As discused in Refs. 
^5),^6|, this model maps onto Anderson's antiferromag- 
netic Ising model with global z quantization axis p8| . 
Consequently, it is normal to recover in Fig. 2 the spec- 
trum of X max (q) found for the pyrochlore antiferromag- 
net |2lJ . Indeed, the largest eigenvalue found in Ref. |2l]] 
is 2 J c ff. Using J c g = J nn + D nn = 5/J/3, we get the flat 
spectrum at A(q)/L> = 10/3, as in Fig. 2. (In Figs. 2-7, 
q is measured in units of 2n/a). The flatness of X max (q) 
reflects the q— space representation of the zero energy ex- 
citation modes in real space. In the kagome Heisenberg 
antifcrromagnet lattice, these are the so-called weather- 
vane modes [§,0. Here, for Ising spins, they correspond 
to closed loops of equal in and out spins which obey the 
ice-rules [|8| ^o). These were first identified by Ander- 
son |2§| ], and have recently been exploited in computer 
simulations of the dipolar spin ice model [42^ . 

In Fig. 3 we see that X max (q) develops an absolute 
maximum for N c = 5. This means that a well defined 
mode becomes soft (massless) at the corresponding in- 
commensurate q or d, which should give rise to incommen- 
surate long range Neel order, unlike what was found in 
Refs. 32, 4J] for their choice N c = 5. We comment on 
this further in Section IV. 

For N c — 10 (Fig. 4), we see that q or d is different 
from its value for N c = 5. Clearly, q or d is not a rapidly 
converging function of N c . We note also that some of 
the flatness present for N c — 1 shown in Fig. 2 is be- 
ing recovered for N c — 10. For N c — 100 (Fig. 5), we 
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clearly note modulations ("ripples") in X max (q) in the 
(111) direction, with now increased flatness restored in 
the overall spectrum. As N c is increased, these modu- 
lations "interfere" in such a manner as to restore a very 
smooth surface, as shown in Fig. 6 for N c — 1000. This 
is similar to what is found in the classical Heisenberg py- 
rochlore lattice with long range dipole-dipole interactions 
where the amplitude of the ripples in the (111) directions 
are continuously decreasing and mutually cancelling out 
each other as the number of nearest neighbor distances 
considered in the dipolar lattice sum is increased to in- 
finity and where, in the limit N c — > oo, a degeneracy 
line occurs at the star of (111) ^2),^). Essentially, the 
ripples get smoothed out and pushed to zero as the finite- 
size effect of the cut-off wavevector q c ~ 2tt/R c (N c ) goes 
to zero. This is more explicitely shown in Fig. 7 where 
X max (hh0) is shown. The convergence at strictly q = is 
slow since the dipolar lattice sum is conditionally conver- 
gent for this q value. This is the origin of the downward 
"spike" in Fig. 7 for q = 0. 




<00i> <hh0> 
FIG. 4. X max (q) vs q in the (hhl) plane for N c = 10. 




<00l> 



<hh0> 



FIG. 2. \ max (q) vs q in the (hhl) plane for iV c = 1. 




<00l> 



<hh0> 



FIG. 3. \ max (q) vs q in the (hhl) plane for iV c = 5. 




<00l> 



<hh0> 



FIG. 5. X max (q) vs q in the (hhl) plane for N c = 100. 




<00l> 



<hh0> 



FIG. 6. X max (q) vs q in the (hhl) plane for N c = 1000 
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<hhO> 



FIG. 7. A maa; (q) vs q for q in the (100) direction, for 
N c — 1 (solid line, no symbols), N c = 2 (right triangles), 
N c — 5 (up triangles), iV c = 10 (down triangles), iV c = 20 
(diamond), N c = 50 (squares), N c = 100 (stars), N c = 200 
(pluses), N c — 500 (crosses), and N c — 1000 (circles). 



IV. DISCUSSION 

Figure 6 constitutes our main result. We see that as 
N c — ► oo, the spectrum of X max (q) becomes very fiat, 
similar but not identical to the nearest neighbor spin 
ice model (Fig. 2). This shows that it is the very 
long distance behavior of dipolar interactions leads to a 
high frustration that restores the symmetry and loops as 
low-energy excitations obeying the ice-rules |^8|,^2| , and, 
therefore, causes the spin ice phenomenon in Ising py- 
rochlore rare-earths systems. One notes that A max (q) 
displays a weak maximum at q or d = (27r/a)[l, 0, 0]. 
The low- lying maximum and weak dispersion of A maa: (q) 
when N c — > oo with respect to the overall smooth and flat 
profile of the spectrum implies a small ordering tempera- 
ture, T c , compared to the nearest neighbor dipolar energy 
scale, D, in agreement with our recent loop Monte Carlo 
simulations where we found T c /D ~ 0.11 Q. Also, q or d 
found here in the mean-field theory for iV c — * oo is the 
same ordering wave vector as the one found in our Monte 
Carlo simulations |fi^| . Consequently, the answer to the 
question raised in the Introduction does appear to be 
that there is indeed a high level of frustration in spin ice 
materials that come from the energetics beyond nearest- 
neighbor. However, this frustration is not perfect, and 
there is indeed a lifting of the degeneracy caused by the 
dipoles, but only at a temperature scale T/D ~ 0.1. 

The results presented here of the quasi-degeneracy in 
X max (q), and the consequential spin ice behavior being 
recovered when N c — > oo (as found in our Monte Carlo 
simulations jll],[l2) also), raise some questions as to the 
validity of the simulation work of Refs. |32| and Hjwhere 
dipole-dipole interactions were cut-off for N c = 5 J32| and 
N c = 12 || . We believe that the results of Refs. §| 
and |44| are invalid as pertaining to real materials due 
to the truncation used for the dipolar sum. Specifically, 
for the values used in Ref. ]32j to model Ho2Ti207, and 
with a choice of N c = 5, we have found in Monte Carlo 
simulations that the paramagnetic scattering is indeed 



incommensurate [ p5[ , in agreement with our mean-field 
calculations. The constraint to work with a small num- 
ber of unit cells and system size L incompatible with an 
incommensurate q or d that would be selected for N c — 5 
for a thermodynamically large sample may precipitate 
the system into a partially ordered state via some type 
of freezing transition as found in Ref. p2fl . Conceptu- 
ally, this is similar to the difficulties in finding the ap- 
propriate long range ordered vortex lattice ground state 
in Monte Carlo simulations of dense frustrated Joseph- 
son junction arrays in a magnetic field, and where glassy 
behavior is found when too small system sizes are consid- 
ered [|||J47|- At any rate, what happens in a simulation 
with N c < O(10 2 ) nearest neighbors is a moot issue in 
terms of modeling and understanding spin ice in real ma- 
terials as our analysis here shows. 

In conclusion, we have shown that one can understand 
at the mean-field level that it is the true long range nature 
of dipolar interactions that cause the quasi-degenerate 
ice-rules obeying states in Ising pyrochlore magnets such 
as Ho2Ti207 and Dy2Ti20y. However, similar to the so- 
called energetic ice models j48|]49|] , dipolar interactions 
do slightly energetically select a preferred state with long 
range Neel order at very low temperature. It would be 
interesting and useful to understand what is the specific 
relationship between the symmetry of long range dipolar 
interactions for local (111) Ising spins coupled via long 
range dipolar interactions, and the symmetry of the py- 
rochlore lattice that leads to an asymptotic restoration 
of a approximately ice-rules obeying quasi-degenerate 
ground state manifold. 
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